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■^ Abstract 

cn 

00 

• Hadronic r decays provide several ways to extract the Cabbibo-Kobashi-Maskawa 

(^ (CKM) matrix element Vus- The most precise determination involves using inclusive r 

CO decays and requires as input the total branching ratio into strange final states. Recent 

'T^ results from B-factories have led to a discrepancy of about S.Aa from the value of Vus 

^ implied by CKM unitarity and direct determination from Kaon semi-leptonic modes. 

K> In this paper we predict the three leading strange r branching ratios, using dispersive 

5-H parameterizations of the hadronic form factors and taking as experimental input the 

measured Kaon decay rates and the r — )■ Kiruj- decay spectrum. We then use our results 

to reevaluate Vus, for which we find \Vus\ = 0.2207 it 0.0027, in better agreement with 

CKM unitarity. 



1 Introduction 

Inclusive hadronic decays of the r lepton provide a unique laboratory to study QCD at low 
energy fl]. However, predicting exclusive decay rates is a notoriously difficult task, that 
requires knowing the relevant non-perturbative form factors over a wide kinematical range. 
While near threshold rigorous chiral perturbation theory (ChPT) methods can be employed, 
the allowed kinematical region extends well into the resonance domain, where different non- 
perturbative tools are needed, such as a combination of dispersion relations and data. 

Focusing on r decays into strange hadrons (see Tab. [l| adapted from Ref. [2]) one notices 
that Tio = Tr-^K-ur^ Tie = Tr-^K-nOur and Tgg = T^-^^-KO^^, which represent 68% of 
the total strange width, are crossed channels from kaon physics. This suggests that, assuming 
lepton universality, one can predict T^-^x-^^, r^-_>.^-^Oj,^ and ^r-^n-ROur using the following 
ingredients: (i) kaon branching ratios (BRs), precisely measured; (ii) shape of the Ktt form 
factors determined by a combined fit to the Ki^ decay distribution and the r~ — )■ KnUr 
invariant mass distribution using a dispersive parametrization for the form factors as presented 
in Refs. [3||4]; (iii) theoretical input on the electromagnetic and isospin breaking corrections. 

The primary purpose of this work is to predict the leading strange r branching ratios 
along the lines outlined above. We will then use the predicted BRs to update the extraction 
of Vus from inclusive r decays ^|6] and explore how this affects the 3.4o" discrepancy with the 
extractions of Vus based on CKM unitarity and kaon decays ff] . 

The paper is organized as follows. In Section |2] we review the prediction of r — )■ Kuj. from 
Kfj^2- In Section |3] we discuss all the ingredients needed to predict r — )■ RttUt branching ratios 
in the Standard Model and give our results and error estimates. In Section |4] we work out 
the implications of the new predicted strange BRs on the inclusive extraction of Vus, and in 
Section |5] we give our conclusions. 

2 T ^- Kiy-T- from K^2 rate in the Standard Model 

Assuming t — n universality in the charged weak current, the r — )■ K^v^ decay rate can be 
predicted from the K — )■ //z/^ decay rate: 

^l Sew f'^-^K/^lV ^r ^r/K 



^Hr^Ku^) = —^^^[- fj-^\ ^ i?^^^^ BR(K,2) , (1) 

2mKml Si^ V 1 - mJ^K J "^K 



fetime respec- 
evaluated 



with Tr = 290.6(1.0) fs {8] and tk = 12.384(15) ns I?) the charged r and kaon h 

t/K — J I 

tively. S'evv represent the short distance electroweak radiative corrections |9||10_ 

at the scale fj, = rrir and /i = mp, respectively. R^m — 1-0090(22) [11] the long-distance 

electromagnetic corrections. Using Eq. ([I]) one finds 

BR(r -^ Kur) = (0.713 ± 0.003) x 10"^ . (2) 



Branching fraction 


HFAG Winter 2012 fit 


Tio = K-Ur 


(0.6955 ± 0.0096) 


10-2 


Tie = K-n^Ur 


(0.4322 ± 0.0149) 


10-2 


T23 = K-2Tx^Ur (ex. K^) 


(0.0630 ± 0.0222) 


10-2 


T2s = K-?>T^^Ur (ex. ir°,r/) 


(0.0419 ±0.0218) 


10-2 


Tss = -n-R^r 


(0.8206 ±0.0182) 


10-2 


140 = 7r-/?°7rV^ 


(0.3649 ± 0.0108) 


10-2 


144 = TX-K^TX^Tl^Vr 


(0.0269 ± 0.0230) 


10-2 


153 = K^h-h-h+Ur 


(0.0222 ± 0.0202) 


10-2 


ri28 = K-r]Vr 


(0.0153 ±0.0008) 


10-2 


ri3o = K-Ti^rfPr 


(0.0048 ±0.0012) 


10-2 


ri32 = -K'K^riPr 


(0.0094 ±0.0015) 


10-2 


Tisi = fC-wz/^ 


(0.0410 ±0.0092) 


10-2 


r8Ol = i^-0Z/.(0^i^K) 


(0.0037 ±0.0014) 


10-2 


r802 = i^-Tr-TT+z/T- (ex. K^^u) 


(0.2923 ± 0.0068) 


10-2 


Tgos = i^-7r-7r+7r°z/^ (ex. K'^,u,r]) 


(0.0411 ±0.0143) 


10-2 


Tiio = X-i/^ 


(2.8746 ± 0.0498) 


10-2 



Table 1: HFAG Winter 2012 Tau branching fractions to strange final states f2j. 



KirVr branching ratios in the Standard Model 



3.1 Relating K -^ irivi and r — >^ KttVt rates 

The decays r — )• KtcUt- and K — )■ ttEui {£ = e, n) are generated by the same underlying 
quark-lepton level operator in the charged current effective Lagrangian (with the replacement 
T ■v^ i). This is true in the Standard Model (SM) and in any extension that respects lepton 
universality. Therefore, the hadronic matrix elements for the above two processes are related 
by crossing. Considering only the SM operator, the K — )■ Triui amplitude involves 



{n{p^)\sjf,u\K{pK)) 






t 



'-{PK-Pn)X^it) + 



{PK + Pn),, 1^^^^ ^ ^"'''^ 



frit) , (3) 



where t = {px — P-k)"^ and ^k-k = ""^i- — f^n- The vector (scalar) form factors /+(t) (/o(^)) 
represent the P-wave (S-wave) projection of the crossed channel matrix element {K7r\s^'^u\0) . 



The scalar form factor /o(t) can be expressed in terms of /+(t) and f-{t) as /o(t) = f+(t) + 
t/AKnf-{t), and by construction, /o(0) = /+(0). The hadronic matrix element relevant for 
r — )■ KfiUr reads 



{K{pK)7r{p^)\s-f^u\0) 



A 



Kn 



{PK+P.)jti 



s - 



{PK - Ptv)^i -{PK + Pn) 



f- 



Kit I 



(4) 

with in this case s = {pK+Pn)"^- The decay rates for r — )■ KnUr and K — )■ Triue involve integrals 
of the form factors over the appropriate phase space. The overall normalization, common to 
both modes is controlled by /_|^'^(0). It is therefore convenient to factor out f^ '^ (0), denoted 
/_i_(0) in the following, in the K^^ and r — )■ Ktip^ decay rates. The phase space integrals 
depend then on the normalized form factors, defined by 



u 



u 



-, Us) 



Us) 



, /+(0) = /o(0) 



With the above definitions for the hadronic form factors, the Ki^ decay rate reads 



T{K ^ 7r£z/47]) 



^^cis^^iKsium) 



'^K 






SU(2) 



y 



(5) 



(6) 



Here 5"^ represents the short distance electroweak radiative corrections m\Wl evaluated at 



the scale /i = nip, Ck the Clebsch-Gordan coefficients, equal to 1 for K^ and l/\/2 for K . The 



12 



13 



quantity 5^-^ encodes the channel dependent long-distance electromagnetic corrections 
and 5^-^(2) ^^® correction for strong isospin breaking. It is defined to parameterize the difference 
between the K ^ -n and K^ — )■ tt^ form factors, so that (5|^(^) = and ^|^(2) 7^ 0. Finally, 
the dimensionless phase space integral is given by 



'X 



dt 



m' 



X"^ 1 + 



K 



2t 



2t 



\Ut)\'+ ^ ^^ 



{2t + ml)\ 



i/owr 



(7) 



with A = [t — {rriK + mT,)'^][t — {rriK — ^tt)^] and tmax = {jriR — ""^tt)^- 

The r — )■ KnVr decay rate has a structure similar to T{K — )■ TidviYA). Including electro- 
magnetic and strong isospin breaking corrections one has 



r(r ^ K-KVr[l]) 



967r3 



/^2 QT 



{\Vus\Um'lK[l + 5^^ + 5f'- 



SU(2)j 



S'eyv represents the short distance electroweak radiative corrections |9|10| evaluated at the scale 
/i = m^. Cx is the Clebsch-Gordan coefficient defined above. 5^{j^ is the channel dependent 
long-distance electromagnetic correction and (5|^^2) ^^^ correction for strong isospin breaking. 
As before, ^|^/^) 



and (5|^(2) 7^ 0. Note that 5^(2) t^ '^su(2) because the K and r decay 



rates involve integrals of the form factors over very different energy regions, 
dimensionless phase space integral , /^ is given by 



Finally, the 



'■K 
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rat 



sr 



ds 



1 



rat 



1 + 



2s 



rat 



<lUs)\h{s)\^ + ^-^ 



gi^^(s)|/o(s)p 



(9) 



with Sk-k = {mK + iTT'ny and qKn the kaon momentum in the rest frame of the hadronic system: 

QKn = — ^a/(s - sxtt) (s - ^a'tt) X 9 (s - sk^) , tx^ = [itik - m^f . (10) 

Taking the ratios of Eqs. ([6]) and (Is]) and multiplying by the ratio of r and K lifetimes, one 
obtains the following relation for BR(r — )■ KnUr) in terms of the crossed channel branching 
fraction BR{K — )■ vrfz/^): 

2 



2m^ 55w /^ (l + ^EM + ^su"(2)) r. 



BR(r ^ iTTTZ/,) = ^^ ^ ^ ^ 3^ BR{K ^ vrez/e) , (11) 



\^-L ^ "em ^ "sU(2)y 



We will use the above formula to predict BR(r — )■ KiiVr)- All the theoretical and experimental 
quantities involving K^^ decays in Eq. ( [Il| are very accurately known W\. The key new 
ingredients are the phase space integrals J^, that require knowledge of the form factors over 
a wide energy range, and the electromagnetic and isospin-breaking corrections relevant to the 
r decays, 5§^ and ^1^2)- In what follows, we describe in detail the evaluation of these three 
input quantities. Before doing that, we make the following general observations about our 
approach: 

• In order to compute /^ (see Eq. (^9 )), we determine /^o'^ i^) ^Y ^ combined fit to the Kij, 
decay distribution and the r~ —^'Kst^'i't invariant mass distribution using a dispersive 
parametrization for the form factors [Slfi]. 



• Calculations of 5^-^ and ^^^(o) 7^ are not as robust as the corresponding quantities for 
K decays, because a rigorous ChPT analysis can only be performed in a corner of r decay 
phase space. However, we will provide in this paper first estimates for these quantities. 
In order to estimate the electromagnetic effects we will use a point-like description of 
pions and kaons, neglecting all structure-dependent effects both in loops with virtual 
photons and Bremsstrahlung amplitudes. For the strong isospin breaking effects, we will 
obtain a rough estimate by using a parameterization of the s dependence of the form 
factor based on a simple resonance model. In both cases we will assign conservative 
uncertainties to the results we obtain. 

One important consequence of the above discussion is that we will be able to predict BR(r~ — )■ 
K^tt'Vt) more accurately than BR(r~ — )■ K'ti^Vt), since the latter involves the poorly known 

"SU(2) • 

3.2 Kit form factors 

3.2.1 Parametrization of the form factors 

To compute the phase space integrals, /^, one needs to know the normalized Kt^ form factors, 
/+(s) and /o(s) in the two energy regions m^ < s < {uik — "m,^)^ (for K^^ decays) and 
{rriK + tUt^Y < s < ml (for r — )■ KnUr)- To this end, a dispersive representation for the form 



factors has been introduced in Ref. |3j. Here we briefly recall the key ingredients of the two 
parametrizations used. For more detailed see Ref. [i]. For the scalar form factor, a dispersion 
relation with three subtractions is written for ln/o(s), one at the Callan-Treiman point and 
the other two at zero. This leads to the following representation for /o(s) 



fo{s) 



exp 






InC + (s - Axn 

Artt 



\nC X' 



nit 



71 



°° ds' 



Ms') 



SK 



(12) 



The two subtraction constants a priori unknown, InC = ln/o((5A'7r) and Aq, the slope of the 
form factor (the third one being fixed since /o(0) = 1, see Eq. (pi)), are determined from a fit 
to the data. 0o('S) represents the phase of the form factor. In the low energy region 5s < Scut 
we use the S'-wave I = 1/2 Kit scattering phase from Ref. 14 . For the high-energy region, 
see discussion below. 

A dispersive representation for the vector form factor /+(s) is built in a similar way [sl 
15-17 . In this case the three subtractions are performed at s = 0. Hence the dispersive 
representation for /+(s) reads: 



A 



exp 



A'-.^ + ^(A'|-A-)^ ' 



ni. 



rat 



+ 



vr 



ds' 



SKt. 



S'3 is' 



le 



(13) 



Use has been made of /+(0) = 1 to fix one subtraction constant. A'_,_ and A" are the two other 
subtractions constants corresponding to the slope and curvature of the form factor. They are 
determined from a fit to the data. As for the phase of the form factor, 0+(s), we parameterize 
it as tan0+(s) = Im/+(s)/Re/+(s) in terms of a model for the form factor /+(s) that includes 



two resonances ii'*(892) and ii'*(1414), with mixing parameter /3, see Refs. [16 19 



his) 



m\* - KK-HKniO) + l3s 



/3s 



D(mK',TK* 



L'(m^.',r^./; 



with 



D{rhR,Tji) 



;^2 



771 



s- KRRe Hx^ris) - i7nRTR{s) . 



In this equation, ttir and Tr are model parameters and Tr{s) and tzr are given by: 



at 



J- R[S — i R^^Y 3 r ~2\ 



kr 



jR 1927tFkF^ 
friR (aKnirhj^))^ 



(14) 



(15) 



(16) 



with aKiris) = 2qKn{s)/y/s. Hk-k{s) is the i^vr loop function in ChPT 
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19 



We emphasize 



here that 7nR and F^ are model parameters and do not correspond to the physical resonance 



masses and width. To find them one has to find the pole of Eq. ( 14 ) or equivalently the zero of 



Eq. ( 15 ) on the second Riemann sheet. Note that this model inspired by the Gounaris-Sakourai 



parametrization l3, 15-22 is built such that the good properties of analyticity, unitarity and 
perturbative QCD are fulfilled. This model is only valid in the r decay region. Therefore we 



will use it for s < Scut ~ '^r- Hence there will be seven parameters to fit from the data: A'^ 
and A'_|_ the slope and the curvature of the form factor and the resonance parameters used to 
model the phase: tjik* and Tk* the mass and decay width of K*(892) and m^,' and F^.' the 
mass and decay width of K*(1414) and /3 the mixing parameter between the two resonances. 
For the high-energy region of the dispersive integrals Eqs. ( 12p!3 ), (s > Scut ~ "^^) the 
phase is unknown and following Refs. l3l[4l[23l|24], we take a conservative interval between 
and 27r centered at the asymptotic value of the phase which is vr. The use of a three time 
subtracted dispersion relation reduces the impact of our ignorance of the phase at relatively 
high energies. The price to pay is that the correct asymptotic behaviour of the two form 
factors is subjected to a set of sum rules derived in [3||4], which is used to constrain our fit 
parameters. 



3.2.2 



Determination of the Ktz form factors from r — >■ Ktzw^ Belle data and K^^ 
analyses 



described in Ref. I3 . We minimize the following quantity: 



X 




We perform a combined fit to the Belle data 25 as well as the Ki^ data [7], along the lines 



jytheo _ jy; 



a2s - oq^ 
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cxp 



+ 
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(y-2v 



InC - lnC^«\ ^ 1 /InC - lnC^« 



A' -A 



iK.. 
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A' - A' 



'K,, 



(17) 



a 



2v 






where N^^^ and a^<:^p are respectively, the experimental number of events and the correspond- 
ing uncertainty in the i^^ bin. The theoretical number of events in a given i bin is 18, 19 



N„ 



theo 



N,,.,h 



1 



tot^ui 



with A^tnt, the total number of events, 6„ 



1 T~^Knu dy/s 

the bin width and TT-^xnu the total decay rate given in 
Eq. ([s]). We fit the first 76 points from threshold skv to sgt ~ 1.51 where our parametrization 
is expected to be reliable. Note that following Refs. 16, 17, 19 we exclude from the fit the 
points 5, 6 and 7 that exhibit a bump which is not present in the preliminary BaBar data [26]. 



We have tested that including these points in the fit amounts to increase the x^ from 60/68 
to 78/71 without any significant changes in the values of the parameters, which remain within 



the error bars. The second term of Eq. (|17|) encodes the constraints coming from K£^ analyses 

We are usins 



where a dispersive parametrization has been used for the form factors 23 



InC 



Kp 



0.2004 ±0.0091, A 



'K,, 



24 



(25.66 ± 0.41) X 10-3 and p(lnC, A', ) = -0.33 from Ref. 



ling 

If 



V represents the covariance matrix. In the minimization we also impose the constraints given 



by the sum rules Eqs. (15) and (18) 



"2^ 



of Ref. 13 
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4 1 
~ ds' 



with a2s 



'>o(s' 



InC 



9,2 ; 



^/2 



's' -A 



Km 



(X2v = A'f — A'i and 



(19) 



^The constraints from the other two sum-rules, see Ref. [310 are not imposed in the fit, see Eq. (17 1, since 
they are automaticaUy satisfied due to the large band taken for (/lo.as and 0+,as- 
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Figure 1: Fit result for the spectrum of t ^ KnUr. The data in black are from Belle Collab- 
The points in green are projected data for a second generation B factory with 
with the same central values of current Belle data and rescal- 
The dashed violet line represents the scalar 



oration '25l 



integrated luminosity of 40 ah 

ing errors according to the expected sensitivity 

form factor contribution. The dot-dashed blue line is the vector form factor contribution and 

the solid red line gives the full result. 
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a 



2v 



IT 



ds - 



-,/3 
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^Ktt 



The results of the fit are presented on Fig. [T] and in Tab. |2] with the correlations between 
the parameters in Tab. Isj Tab. [2] display results for the fit to real data f25l and also projected 
data from a super-B factory, obtained by keeping the same central values of current Belle 



data 25 and rescaling the errors according to the expected sensitivity of a second generation 
B factory assuming an integrated luminosity of 40 ab~^, see e.g. Re f |27| . Using these results 
we can compute the phase space integrals Eqs. (ItI [qI) given in Tabs. |4| and [SJ 



3.3 Electromagnetic effects in r — >^ KttUt- 

While the electromagnetic corrections are known for Ki^ to order (e^p^) in ChPT |12|[l3|[28| , 
they have never been computed in the case of r — )■ KirUr. In this case there are no rigorous 
methods to compute electromagnetic effects over the entire phase space, because the kinematics 
of r decays allows the hadronic invariant mass squared s = {pK+PnY to extend well beyond the 
chiral regime, all the way to s = m^. While resonance-model calculations are possible 29,30 





T —)■ KnUr k Ki3 


r -)■ KnUr k K^ 




Belle 


2nd generation B factory 
(projected) 


InC 


0.20352 ± 0.00890 


0.19880 ±0.00498 


X'o X 103 


13.824 ±0.824 


13.703 ±0.521 


rfiK* [MeV] 


943.59 ±0.58 


943.76 ± 0.06 


Tk' [MeV] 


67.064 ± 0.846 


67.290 ± 0.088 


m^*' [MeV] 


1392.2 ±57.6 


1361.7 ±6.3 


f;^.. [MeV] 


296.67 ±160.28 


254.62 ± 17.45 


/3 


-0.0404 ± 0.0206 


-0.0338 ± 0.0023 


AV X 10^ 


25.621 ±0.405 


25.601 ±0.277 


XI X 10^ 


1.2221 ±0.0183 


1.2150 ±0.0090 


xV^-o-/ 


60.2/68 


28.1/71 



Table 2: Results for the Kit form factors parameters from a combined fit to t ^^ K-kUt 
and Ki3 . Note that ffiR and Tr are model parameters and do not correspond to the physical 
resonance masses and width. 



Parameter 
InC 

A'o 
rfiK* 



A' 



\nC 


A'o 


niK* 


Tk' 


m^.' 


r^.' 


/3 


K 


K 


1 


0.943 


-0.093 


-0.117 


0.047 


0.005 


-0.003 


0.342 


0.135 


- 


1 


-0.066 


-0.068 


0.040 


0.027 


-0.067 


0.318 


0.266 


- 


- 


1 


0.951 


0.196 


0.240 


-0.345 


0.001 


-0.250 


- 


- 


- 


1 


0.145 


0.179 


-0.273 


0.017 


-0.160 


- 


- 


- 


- 


1 


0.926 


-0.842 


0.088 


0.030 


- 


- 


- 


- 


- 


1 


-0.917 


0.088 


0.030 


- 


- 


- 


- 


- 


- 


1 


-0.128 


-0.018 


— 


— 


— 


— 


— 


— 


— 


1 


0.735 



Table 3: Correlations between the parameters of the fit. 



here we will give a first estimate of the long-distance electromagnetic corrections to r — ?■ KtcUt- 
based on point-like mesons and leading Low bremsstrahlung contributions, i.e. neglecting 
structure dependent effects. With these approximations we provide the corrections to both 
differential and total rate for the processes r — )■ KttVt-. 

The leading 0{a) long-distance EM corrections arise from one-loop corrections to the de- 
cay amplitudes and real photon emission. Only the one-photon-inclusive decay rate is infrared 
(IR) finite to 0{a). Our approach here relies on the analysis of EM corrections io K ^ -nivn 



and r — )■ ttttz/t- presented in Refs. 12,13 and 29,30 , respectively. Adapting the arguments 



Integral 


result 


error 


exp 


tlieo 




0.50418 


0.01762 


0.01689 


0.00501 




0.15472 


0.00022 


0.00022 


0.00000 


^Ro/^RO 


3.25864 


0.11115 


0.10634 


0.03235 


Ir^ 


0.52387 


0.01958 


0.01889 


0.00515 


Ir^ 


0.15909 


0.00025 


0.00025 


0.00000 


^R+/^R+ 


3.29282 


0.12032 


0.11589 


0.03235 



Table 4: Phase space integrals for the charged and neutral modes of r ^ Knu and K^s as well 
as their ratio using the results of the fits to Belle and Ki^ data, see Tab. [^ The experimental 
uncertainty comes from the uncertainties from the fit parameters and the theoretical uncertainty 
comes from the uncertainty of the phase of the form factors in the inelastic region, where a 
large hand of 2tt has been taken, see section 2.3.1. The two uncertainties have been summed 
in quadrature to give the final one. 



Integral 


result 


error 


exp 


theo 




0.49590 


0.00820 


0.00662 


0.00484 


^RO 


0.15471 


0.00015 


0.00015 


0.00000 


^Ro/^KO 


3.20545 


0.05060 


0.03562 


0.03130 


Ir^ 


0.51536 


0.00858 


0.00631 


0.00498 


Ir^ 


0.15908 


0.00017 


0.00017 


0.00000 


^R+/^R+ 


3.23973 


0.05114 


0.03635 


0.03132 



Table 5: Phase space integrals for the charged and neutral modes of r ^ Knu and K^s as well 
as their ratio using the results of the fits to the projected 2^'^ generation of B-factories and 
K£3 data, see Tab. [^ 



presented in Ref. 12, 13 we find that long distance EM effects in r — > KtcUt induce [j 



(i) An overall correction grad{s,u) to the differential decay rate, that combines the effect of 
soft real photon emission and the universal soft part of one-loop diagrams. The virtual- and 
real-photon corrections are IR divergent and depend on the IR regulator M^, while their sum 
is finite: 

rv _ _ _ _ _ _ 

(21) 



grad[S,U) 



a 

2^ 



Tc{u, ml, ml, M?) + ghremsis, u, ml, ml, M^) 



^For the two decay modes we adopt this conventions for the particle four-momenta: t~ (p-r) — ^ 
tt~{pi)K°{P2)i^t{q) and t~{pt) — >■ K~ {pi)!^^ {p2)vt{<i) ■ The EM corrections involve the Mandelstam vari- 
able u — (pt — pi)^, where Pt and pi denote the four-momentum of the r and the charged meson (K or tt) in 
the final state. Moreover, ml = p\ denotes the mass squared of the charged meson. 



The expression for rc(u, m^, m^, M^) can be found in Ref. 
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pleteness in Appendix pi gy,rems{s,u,ml,ml, M^) encodes theBremsstrahlung effects in the 



13 



and is reported for corn- 



leading Low approximation and its expression can be found in Ref. 29 , 30 and Appendix C 



(ii) Shifts to the form factors: f±Q{s) — )■ f±Q{s) + Sf^Q^s, u). These shifts arise already when 
treating K and vr as point-like as soon as one uses momentum-dependent vertices for the weak 
hadronic current. Sf±{u) are given by 



5/ 



K- 



(U 






{u) 



a 


1 


A-K 


/+(0) 


a 


1 



[ri(M, ml, m]^) ± T2{u, ml, m\)\ + . . . , 
\r2{u,ml,ml)±Vi{u,ml,ml)\+ ... , 



(22) 



(23) 



4vr/+(0) 

The dots denote structure-dependent corrections that are hard to estimate over all the phase 
space. Near threshold, the ChPT expressions in terms of low-energy constants can be found 



in Ref. p^ 13 . The loop functions Fi 2(M,m^,'nif) can be found in Ref. 12 13 and in Ap- 
pendix p3l Finally, in terms of the shifts Sf^'^{u), the corrections to the scalar form factor 



cKtti 



U . 



reads 6f^^s,u) = 6f^^u) + s/Ak. 6f_ 

With the above prescriptions, and linearizing in the corrections to the form factors, we 
obtain the following expression for the photon-inclusive double differential rate r — > KtiVt-Yj] 
decay: 



d^T^K-Ku[-t] _ GpC]^ S^^;^\jj^{^)Vus 



ds du 1287r'^m^ 



D^^s, u) (\U{s)\^ + 2Re [U{s)5fl{u)\ ) 



u 



(24) 



+ D^^{s,u) Re 



/+(^)/o (^) + h{s)6f;{s, u) + 5h{u)f*{s) 



X 



5'rad(S,UJ 



The expression for the Dalitz plot kinematic densities D+_o,+o('5,m) can be found in Ap- 
pendix |A} Integrating over the u variable we obtain the EM-corrected distribution in the 
Ktt invariant mass: 



dV KT^Yy] _G],C\S'Ew\f+{^) Vus\ m 



2^3 



ds 



QGtt^Sa/s 



m: 



1 + 



2s 



m.t 



qlAs)\U{s)? [I + ^Uis)] 



As 



QkM \Us)\' [I + 5°m(^)] + Re [h{s)n{s)\ 5ll{s 



(25) 



with 



^em(s) 



!uZS ^^ ^+("'^) {\h{s)\^ 9r^{s,u) + 2Re [/+(5)5/;(w)]) 



C:^ du D^{s,u)\h{sW 



(26) 
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Figure 2: Correction factors (^em(s) (left panel) and S^ui^) (^ight panel) to the vector and 
scalar contribution to the differential decay rates of both r — )■ KnUr modes. 



SUs) ^ 



j;--;;^ du Do{s,u) {\fo{s)\'grUs,u) + 2Re [fo{s)6fS{s,u)]) 



/:it) du D,is,u)\foisW 



3S\/S /'"m='''('') / _ _ 



■'EM^ 



Ami 



"^mini^J 



u 



+ Re[U{s)5f;{s,u) + 5U{u)f;{s)\ 



(27) 



(2^ 



■Umin.maxls) caii be found in the appendix. The functions 5em('5) ^.re shown on Fig. [2| Further 
integrating over the s the distribution (25) with and without electromagnetic corrections, and 



taking the ratio, we get 5^{/[- Assigning an uncertainty of ~ a/vr to the unknown structure- 
dependent corrections, we get: 



rf-7 



-(0.2 ±0.2) 



5- 



ROt 

EM 



-(0.15 ±0.2)% 



(29) 



Note that a comparison between the leading Low approximation and the full calculation is 



performed in Ref . 29 , 30 , and it shows that it leads to a comparable correction to the decay 



rate. Hence we expect this calculation to give a reasonable estimate for the electromagnetic 
corrections to the r — )■ KttVj. total decay rates. We have introduced the EM correction factors 
'^em'*' ("^) i^ ^^^ fitting procedure and we have found that these corrections do not affect the 
determination of the form factors at the current level of precision. 



3.4 Isospin breaking corrections in r — >^ i^ it v-r 

In order to estimate strong isospin breaking effects, we focus on the dominant vector form 



factor. We adopt a simple parameterization of the ratio f^ '^ is)/f+ '" (s) based on a single 
vector meson resonance exchange. The ratio f^l '" {s)/f^ '^ (s) differs from unity because of 
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(i) 7r° — ?7 mixing and (ii) possible isospin-breaking effects in the coupling of K* to Kt:. To 
leading order in isospin-breaking, the first effect is independent of s, and completely controlled 
by the 7r° - 77 mixing angle e = ^ ^^_"2(X+md) = 1-16(13) x 10"^ 31 . The second effect can 
be estimated by using couplings of vector mesons to Goldstone Bosons that involve insertions 



of quark mass matrices, such as those introduced in Ref. 32 . Requiring that the form factors 
in the isospin-symmetric limit fall off as 1/s, single vector meson resonance exchange implies 
the parameterization: 

/r-°W//f'-W-(l + V3.)(l.«i£|j,^.). (30) 

The only unknown parameter in the above expression is the coupling g ~ 0(1), which we 
vary between —2 and +2. This gives a first rough estimate of the effect of s-dependent isospin 
breaking effect, namely 5|^(2) ~ ±0.5%. On the other hand, the constant part due to 7r° — 77 
mixing is better known and is 100% correlated with the analogous i^^s quantity. Putting 
the two ingredients together, this procedure leads to S^-^,^) = (2-9 ± 0.4mixing ± 0.5)%. We 
emphasize that this is a far-from-complete estimate of strong isospin breaking effects, and it 
is only meant to provide a rough estimate of the central value and uncertainty associated with 
these effects. 

3.5 Branching ratios 



Using Eq. (11) we predict Br(r -^ K r^^Ur) and Br(r — t- K^n z/,-). In Tab.p^we summarize 



the input values used for the predictions. We find for the branching ratios 

BR(r" -^ K^-K-Ur) = (0.8566 ± 0.0299) ■ 10"^ (31) 

BR(r" -^ K-7i%r) = (0.4707 ± 0.0181) • 10"^ (32) 

with a 100% correlation. The error comes exclusively from the uncertainty on the r phase 
space integrals. In Tab. [T] results for the 2^^*^ generation of flavour factory with the error 
budget can be found. One sees that the uncertainty coming from the evaluation of the phase 
space integrals can be reduced by a factor of three. Then the uncertainties coming from EM 
corrections start to matter. 



4 Implications for the inclusive determination of V^ 



us 



The most precise determination of \Vus\ from r decays comes from the measurements of inclu- 
sive IAS"! = and [AS*! = 1 tau decay widths. Indeed one can build the theoretical quantity 



Rt,NS Rt,1 

^ud\ 



5Rr,th - TTT-r^ - TTT^ , (33) 



where Rr is defined as 



r[r -^ hadrons i/^] 
Rt = — Wf — ] • l34j 
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Parameter 



Value 



Ref. 



Parameter 



Value 



Ref. 



BR(i^el) 



0.05078(31) 

^12.384 ± 0.015)ns 

(290.6 ± 1.0) fs 



BR(K 



Le3) 



0.4056(9) 
(51.16 ± 0.21)ns 
(290.6 ± 1.0) fs 



CK 

Jew 
S 



EW 



<5.^^ 



(%) 



EM 

"jiV{2) 

''SU(2) 
JT I je 



1.0232 ± 0.0003 
1.0201 ± 0.0003 
0.050 ± 0.125 
-(0.2 ± 0.2) % 
0.029 ± 0.004 
0.029 ± 0.006 



I 

10 



cK 

'-'EW 

S 



EW 



28 
Sec. 





(%) 



3.3 



Sec. 



3.4 



10.32059 ± 0.48240 Tab. 



"EM 

'^EM 
i;K-K 

s.K-n 
''SU{2) 



1.0232 ± 0.0003 

1.0201 ± 0.0003 

0.495 ± 0.110 

-(0.15 ± 0.2) % 





10.21432 ± 0.43058 



10 



28 
Sec. 

0[ 

Sec. 



3.3 



33 



3.4 



Tab. 4 



Table 6: Input used to compute r -^ K ix^Vj. and r — ?■ K^n v^ branching ratios. 



Mode BR % err 

T- -^ K^Ti-Vr 0.8427 ± 0.0122 1.45 
T- -^ K-n^Ur 0.4631 ± 0.0079 1.71 



0.22 
0.06 



tk 
0.41 
0.12 



0.34 
0.34 



Jk/3 

1.24 
1.25 



EM SU(2) 
0.46 
0.47 1.00 



Table 7: Prediction for the r^ — > K^n^Uj. branching fractions in % from the Kf.3 branching 
ratio using the 2^'^ generation of B factory projected results for the phase space integrals, see 
Tab. \^ The different sources of uncertainty are given. They have been summed in quadrature 
to give the total one. 



This quantity vanishes in the 5'f/(3) limit and can be precisely determined within QCD 
combining perturbative QCD and low energy data (5l|6,34 



Hence, we can extract Vus 
from Eq. (33) using the theoretical estimate of SRr,th and the precise measurements of non- 
strange {Rt,ns) and strange {Rt,s) inclusive decays, and \Vud\- Following Ref. [2], we take 
SRr^th = 0.240 ± 0.032, with a systematic error on \Vus\ that lies between the two more recent 
estimates |35|[36|. We use |Kd| = 0.97425 ±0.00022 from the superallowed 0+ -^ 0+ nuclear /3 
decays 37 . Using the HFAG Early 2012 averages from the r branching fractions reported in 
Table [l together with their reported statistical correlations and replacing the r — )■ Kur and 
r — 7- KirUr results by our prediction in Eq. ^ and Eqs. ( [3l] ), (32), we obtain [2] 



BRr^s = BR(r -^ X^Ur) = (2.9648 ± 0.0656) • 10 



-2 



(35) 



With this estimate and BR^,^^ = (61.85 ± 0.11)% and BR^,e = (17.839 ± 0.028)% |2j, we get 

iKs I = 0.2207 ±0.0027. (36) 
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Figure 3: Determination of \Vus\ from semileptonic, leptonic kaon decays ^, hyperon de- 
cays ^381 and inclusive and exclusive t decays ^. The errors bars correspond to the determi- 
nation from exclusive r decays (blue), the inclusive hadronic r decays (red), and our prediction 
(cyan). The grey band displays the value of \Vus\ assuming unitarity of the first row of the 
CKM matrix. 



We summarize in Fig. p^the different extractions of \Vus\ from semileptonic and leptonic kaon 
decays, hyperon and r decays. Our prediction shifts the inclusive determination of \Vus\ 
towards the exclusive one by ~ Icr. 



5 Conclusion 

The experimental precision of data on leptonic and semileptonic kaon decays matched by sub- 
percent theoretical calculations allowed the most accurate determination of Vus ||7|. Assuming 
lepton universality, we use the same data in combination with the measurement of the K% 
invariant mass distribution in the r — )■ Kttv decay and a dispersive parametrization for the 
form factors [3||4] to obtain a precise prediction for about 68% of the total strange width. A 
first evaluation of electromagnetic and SU(2) breaking effects has been derived to this purpose. 
We find: 

BR(r- ^ K-Pr) = (0.713 ± 0.003) x 10"^ , 

BR(r- -^ K\-i^^) = (0.8566 ± 0.0299) ■ lO'^ , 

BR{t- -^ K-n^u^) = (0.4707 ± 0.0181) ■ 10"^ , 
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^3 = BR(r -^ Kur) + BR(r- -^ K^tt-u^) + BR{t- -^ R-n^Ur) = (2.040 ± 0.048) ■ 10"^ . 

B^ is 1.7(7 higher with respect to the world average measured value. In addition we obtain a 
determination of Vus from inclusive tau decays using the above prediction for the branching 
ratios and the world average values for the rest of tau branching fractions. We find: 

\Vus\ = 0.2207 ±0.0027, 

and for the unitarity of the CKM quark mixing matrix as applied to the first row, we obtain: 

1 - iKsT - iK^r = 0.0021 ± 0.0013 (-1.6a) . 

Finally, we have shown that measurements of the Kn invariant mass distribution at a second 
generation B factory with integrated luminosity of 40 ab~^ would reduce the uncertainty in 
the T — )■ KttUt- BRs by a factor of three, and therefore further reduce the error on V^s- 

A Kinematic densities 



The differential decay rate Eq. (24) involves the kinematic functions 



2 

D^''{s, u) = —^{ml - s) + 2m\ml - 2u {ml - s + ml + ml) + 2u^ 

-ml{2u + s-ml-2ml) + -^'^{ml-s) , (37) 

s s'^ 2 

fl-,.,.).^(l_^), ,38) 

D^;(s,u)=^i^ Uu + s - ml - 2ml - ^ « - s)\ , (39) 

with A21 = m^ — m^l. The above expressions are valid for both decay modes, with the follow- 
ing conventions for the particle four-momenta: T~{pr) — )■ 71~{pi)K^{p2)i't{<i) and T~{pr) — )■ 
-^~(Pi)^°(P2)i^r('?)- The Mandelstam variable u = (Pr—pi)^, where Pr and pi denote the four- 
momentum of the T and the charged meson {K or tt) in the final state. Moreover, mf = pi 
denotes the mass squared of the charged meson. 

B Loop functions 

We now give expressions for the loop functions characterizing virtual photon corrections. We 
denoted by M the charged meson mass, by m,i — )■ rrir the charged lepton mass, and by M^ 
the photon mass used as infrared regulator. In order to express the loop functions Ti^2,c in a 
compact way, it is useful to define the following intermediate variables: 



m; _ ^ „ V ,, Y - VY^ - AR 



,2 
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In terms of such variables, of the dilogarithm 



Li2(x) 



^ dt 



log(l-xt), 



(41) 



and the auxihary functions 

1 X 



2 i\/r2\ 



C{v,mj,M 



n' 1 



+ Lt2{X^) + Lt2 ( 1 - ^ ) + Lt^il - xVR) 



J^{v,mj,M^) = -j= . ^„„ InX 



Rl-X^ 



we have: 



o o Ox o , o Ox Mmp ( Xy loeX 



(42) 
(43) 

(44) 



and 



ri(t;,m,^M^) 



Y^^v.mlM^) 



- Ini? + (4-3r)J^(t;,m2,M2) 

7T7 \ r "I I 

1-^) -^(t;,m2,M2)(l-i?)+lni? - -(3 - F)^(t;,m2, M^) .(45) 



1 /^ m^ 
2 



C Real photon emission 



In Refs. 



29 



30 



the function (y'brems(s, M, "n^i, '?«2) ^7) is denoted by 5'brems('S, m, M^), omitting the 
dependence on the meson masses vn?^^. We report here the full expressions for completeness. 
The Bremsstrahlung function is given by: 



a 



ghrems{s,U,M^) = -[Ju{s,U, M^) + J2o{s,U, M^) + Jo2{s,U, M^)] 



TT 



Jn(.,.,M,) = log(^i^j^log^^-^ 



(46) 
(47) 



+ ^ {Li2{l/Y2) - Li2{Y^) + log2(-l/r2)/4 - log2(-l/ri)/4} (48) 



J2o(s,M, M^) = lo^ 



Jo2{s,u,M^) = loj 



M^{ml - s) 

ITIt-X+^S, u) 



M^{m1 + m\ — s — u) 
mix^{s, u) 



(49) 
(50) 
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The auxiliary variables are: 



x±{s,u) 



2m{ 



2 ml {ml + s) -{s + ml- Mj) {ml + ml - u) 



± v/A(s,m^,m2) A(u, mf, m^) 



Fi 



1,2 



a 



P 



7 



1 - 2a ± v/(l - 2a)2 - (1 - /32) 



(m^ — s)(m2 + ml — s — u) A(u, ml, m,l) 



{ml + m^ — u) 
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^yX{u,ml,ml) 
m^l + m^ — ^i 

^/X{u,ml,ml) 
2VI 



5 = —m^^m^l + mbl{ml — s){mb\ — u) — su{—m,l + s + u) 
+ml{—m.^ + SU + r 
with \{x, y,z) = x^ + 1/2 _|_ ^2 _ 2 (^j^y _|_ 3,^ _|_ ^^^^ 



+ml{—m.^ + SU + m^s + m^w) 



(51) 
(52) 

(53) 
(54) 
(55) 



(56) 
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